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A Reply to Mr. Nelson's Strictures on Delaunay's Method of De¬ 
termining the Planetary Perturbations of the Moon. By 
G. W. Hill. 

For several years past Mr. Neison has been maintaining in 
the Monthly Notices and Memoirs of the Society that Delaunay’s 
investigation of the two long period inequalities in the Moon’s 
motion arising from the action of Venus is seriously defective, 
on account of the omission by him of a certain class of terms. 
In the Monthly Notices for last June there appears a long article 
by him upholding this view; to this I wash more especially to 
direct attention. 

At the outset I may be allowed to say that all this criticism 
is without foundation. It appears to arise, partly from the very 
confused conception Mr. Neison seems to have of the nature of 
Delaunay’s method, and partly because he fails to notice that 
Delaunay, after setting the degree of approximation he wishes 
to attain, always rigorously adheres to it. If we were obliged 
to admit the validity of all the statements in this article, an easy 
oorollary from them would be that Lagrange’s general method 
of the variation of arbitrary constants in the problems of 
mechanics was a blunder. Now, I think that no one acquainted 
with this method could, for a moment even, entertain such a 
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2 Mr. Hill, Reply to. Mr. Nelson's xlyii. i, 

proposition. Hence we may conclude there is some flaw in the 
reasoning of this Paper. But this must be substantiated by 
noticing seriatim the objectionable points. 

In the first place, why bring forward Hansen’s published 
values of the coefficients of these inequalities for the purpose of 
throwing discredit upon Delaunay’s values, when their author, 
himself, virtually confesses, be has no confidence in them, by 
saying he had computed them in two different wavs, and found 
essentially different results F And, to the very end of his life, 
he appears never to have been able to find out whether one of 
these results was right, and which it was, or whether both were 
wrong. It would be an amusing circumstance should it turn 
out that the set of values, withheld from publication by Hansen, 
were identical with those of Delaunay. 

There is some inexactitude in Mr. Heison’s statement, re¬ 
garding the degree of approximation adopted by Delaunay in 
calculating the coefficient whose argument is SI" — i$V. In this 
connection we note that, on account of the close proximity of 
the Moon to the Earth, a planet cannot produce in her motion 
inequalities of the same order with those it produces in the Earth, 
but that they are only of the order of these multiplied by the 
solar disturbing force ; and this is as true of the indirect action 
as of the direct. How, on referring to Delaunay’s work, we see 
that he has considered, not only the term of the lowest order in 
each portion of the coefficient, but also multiples of this by m 2 . 
Hence, it is correct to say that he has considered terms of the 
order of the mass of Venus multiplied by the square of the solar 
disturbing force, but not those multiplied by the cube of the 
latter. 

Mr. Heison regards the evidence adduced in his earliest Paper, 
as conclusively establishing the omission by Delaunay of a certain 
class of terms. But what was this evidence F Simply, that 
Hansen was at variance with Delaunay. How, since Hansen 
was as much at variance with himself as he was with Delaunay, 
what weight ought to be attributed to this evidence F Then, 
Mr. Heison believes that certain discrepancies between results, 
obtained on the one hand by himself, and, on the other, by M. 
Gogou and myself, have their origin in the same cause. How, 
if Hansen’s investigation and also Mr. Heison’s were accessible, 
this point could be immediately pronounced upon; but since 
they are not, it appears useless to speculate on the matter. 

Mr. Heison says (p. 416), “He (Delaunay) substitutes the 
preceding value for the term in the disturbing function with the 
argument 4 in the differential equation and integrates.” This 
does not correctly represent what Delaunay does. For he sub¬ 
stitutes in the differential equation not only the term factored by 
cos Z, but also the non-periodic portion of It: to wit, in the 
memoir of 1862, the terms 


A. 

2 a 




a 2 ri+i 

a' 3 L 4 8 


c* + 


225 

64 


n V32 32 )ri l J 


© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 


Downloaded from http://mnras.oxfordjournals.org/ at Florida International University on June 10, 2015 




1886MNRAS..47....1H 


Nov. 1886. 


3 


Strictures on Delaunay's Method. 


and, in the memoir o£ 1863, the terms 
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The terms differ in the two cases, because the degree of 
approximation aimed at requires the preservation of different 
terms with allowable neglect of all the rest. From this non¬ 
periodic portion of ft results, in both cases, much the larger part 
of the two inequalities considered. Mr. Neison’s failure to note 
this completely invalidates his argument on the two following 
pages, by which he attempts to prove the incompleteness of 
Delaunay’s procedure. And, in this connection, it may be noted 
that it is not necessary that the coefficient in (38) should vanish 
identically in order to prove Delaunay right; it is necessary 
only that it should turn out of such an order of smallness as to 
prove that the adopted degree of approximation had been 
attained. 

On p.417 it is said that the coefficients B and B' “only 
differ by small quantities, unimportant for the present purpose.” 
So far from this being the case, the difference B' —B constitutes 
one portion of the terms which Mr. Neison, all along, has been 
asserting were neglected by Delaunay. 

That Delaunay, in treating the two Venus inequalities, dis¬ 
carded his own method, and employed the old one recommended 
by Poisson, is erroneously stated on p. 423. The fact is that 
the method followed is the same as that he had used in deriving 
the solar perturbations. Next, Delaunay is found fault with 
(p. 424) because he confines himself to calculating in ft the term 
which has the argument of the particular inequality he is dealing 
with; while it is plain that there are a multitude of terms in ft, 
having other arguments, which could contribute to the value of 
the coefficient sought. This is true, but Delaunay’s reasons for 
passing by these terms are quite evident, fn the first place, it 
must be remembered that his final expression for the inequality 
is a formula of substitution, which must be made, not only in 
the mean longitude of the Moon, but also in the equation of the 
centre, in the eveetion, variation, and in all the inequalities 
arising from solar action. Hence, Delaunay’s method of treat¬ 
ment enables him to obtain, with very little additional labour, 
all the terms in the expression for the true longitude which 
involve the very small divisor arising from the slow motion of 
the argument which he is considering; and that whatever may 
be their arguments. And, secondly, while the terms in R, 
having other arguments, which would be treated by Delaunay 
as giving rise each to a distinct transformation, can, in a strict 
sense, add something to the coefficient of the inequality in the 
true longitude, practically these terms are insensible ; for, al¬ 
though they may be of the same order, before integration, as 
the quantities retained, they are altogether independent of the 
excessivelv small divisor which arises from the slow motion of 
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the argument of the inequality. As illustrating this point, it 
may be remarked that, in the case of the two Venus inequalities 
in question, we get such relatively large coefficients as 16" and 
o //# 2 7 only by multiplying the corresponding terms in R by 
factors which are about 15,000,000 in the first, and 10,000,000 
in the second inequality. Hence, if there are other terms, which 
rigorously ought to be added to the preceding values, but which, 
while in other respects of the same order of smallness, have 
factors not much exceeding unity, it is very apparent they may 
be neglected. 

In the next place we find Delaunay charged with neglecting 
every term of the solar perturbations save the term of the lowest 
order in the variation in calculating the proper form for R. 
And it is said that his development “in no sense depends on his 
method of transformed elements, though made to appear as if it 
does; nor does it differ in anyway from the values hitherto 
employed by astronomers save in being somewhat less complete.” 
These statements misrepresent Delaunay. He arranges under 
four different heads the transformations made by him, and they 
involve no less than 16 out of the 57 operations of his first 
volume, besides 4 complementary ones. And whether the amount 
of work in this be regarded as much or little, I have ascertained 
that it is precisely sufficient to obtain the degree of approximation 
he proposes in the coefficients B, viz., to terms involving m 2 . 
Carrying the approximation farther could only have afforded 
him terms of a higher order. It is, of course, open to Mr. 
Heison to say he deems this degree of approximation insufficient; 
and nothing can be said in opposition. But this is very different 
from saying Delaunay has committed errors. Again, I am not 
aware of the existence of any published investigation in which 
the degree of approximation is greater. 

The reasoning Mr. Heison employs to show that Delaunay 
deserts, in this investigation, his own method and returns to the 
old method recommended by Poisson, is certainly very strange. 
He notes that the differential equation used has nothing in it to 
distinguish it from the corresponding one which Poisson would 
have used. But from what circumstance does this state of things 
arise P Simply because it is Delaunay’s habit to omit, in the 
statement of his equations, every term which gives rise, in the 
final result, only to terms of a higher order than he has agreed 
to retain. The factors in question, in Delaunay’s method, can 
be expressed only as infinite series; it is necessary, therefore, to 
cut them off at some point, and he determines this point in the 
way just stated. If reference is made to the same equation, in 
the memoir where Delaunay treats the other Venus inequality, 
it will be found to be duly distinguished by the presence of 
additional terms, Delaunay writing as many as are just sufficient 
for his purpose. 

Mr. Heison next notices two assumptions, which he says have 
been made by Delaunay in his integration. 
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The first is that the factor 2 , which multiplies —, is treated 

an di 

as if it were constant. But here he forgets that, with. Delaunay, 
at this stage of the work, the symbols a, e. y, 7 , g, and h, denote 
quantities which have no solar perturbations; and that, eonse- 

quently, the deviation of — from a constant has the mass of the 

an 

planet as a factor. Thus, as already has this factor, the 

additional terms, which would in this manner arise, would have 
the square of the mass of the planet as factor ; these, as all ether 
investigators, Delaunay expressly neglects. 

With regard to the second assumption, in reference to which 
Mr. Keison makes what he thinks his chief point against 
Delaunay, let us consider what is the essential difference between 
Delaunay’s method and that employed by the earlier investi¬ 
gators. Delaunay said to himself, Do not let us go back to the 
elements of the Keplerian ellipse every time we have to consider 
the action of a new force on the Moon, but let us determine our 
new wave of motion in such a way that it may be superposed on 
the curve which the Moon would describe under the action of all 
the forces pi’eviously considered, instead of on the Keplerian 
ellipse. At any stage of progress, in expressing the Moon’s 
co-ordinates, there must, of necessity, appear in them six arbitrary 
constants which have been introduced by integration. Let us 
take these as variables, instead of the six elements of the Kep¬ 
lerian ellipse. This course demands that the differential equations 
employed by the earlier investigators should be somewhat modi¬ 
fied. The modification appears as a change in the values of the 
quantities which Poisson denoted generally by the symbol [<x, &]. 
Kow, just as it would be absurd to maintain that the elements 
of the Keplerian ellipse suffer perturbations from the action of a 
centrobaric Earth, so it is absurd to maintain that the quantities 
a , e , y, 7 , g and h, employed by Delaunay after he has got through 
with the solar perturbations and has arrived at the treatment of 
the planetary perturbations, and which are the elements of the 
curve which would be described by the Moon under the combined 
action of the Earth and Sun, suffer perturbations from the latter 
body. Yet Mr. Neison’s argument, when divested of its ob¬ 
scurities, is seen to be nothing more or less than a. plea that 
these quantities do suffer perturbations from the Sun. 

To make the matter plainer, let us suppose that Delaunay, 
groping about in the dark, had fallen upon the Poissonian 
equations, and, thinking them to be his own, had used them as 
such; and, moreover, on making his substitutions, had made 
them only in the elliptic portion of the co-ordinates. Then he 
would have committed the very error Mr. Xeison lays to his 
charge. But since he uses equations suitably modified to the 
new signification of the quantities a t e, Ac., and, moreover, 
makes his substitutions in the complete expressions for the 
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Moon’s co-ordinates, and not in the elliptic portion only, as the 
earlier investigators do, is it not plain that, by these two modi¬ 
fications, he obtains terms which he would not have obtained in 
the former supposed case ? Now these terms, in sum, are pre¬ 
cisely equivalent to those Mr. Neison accuses him of neglecting 
by omitting to include B/" in his disturbing function. Thus it 
is seen that Delaunay takes account of ~SL r,! in an indirect manner, 
the peculiar nature of his method absolving him from considering 
the terms arising from B/" as a separate class. 

Perhaps the matter will be clearer still if we say that, just 
as in determining the solar perturbations we have no class of 
terms of the order of the product of the mass of the Earth by 
the mass of the Sun, simply because the Earth’s action is con¬ 
sidered as the principal force, so when we come to treat the 
planetary perturbations by Delaunay’s method, there is no 
special class of terms of the order of the product of the Sun’s 
mass by the planet’s mass, for the reason that here the combined 
actions of the Earth and Sun are regarded as forming the prin¬ 
cipal force. 

Next we must not pass over without notice the quite 
erroneous method Mr. Nelson proposes (pp. 430, 431) for getting 
the proper expressions for the Poissonian quantities [a, &] ; 
viz., by substituting for the elements in the expressions proper 
to the older form of the differential equations their complete 
values as functions of the time, and then neglecting all the 
periodic terms. It is very certain this procedure will not give 
the same values as Delaunay has, who obtains them by taking 
the partial derivatives of a , e, and y with respect to the elements 
L , G, and II, which are the conjugates of l, g , and h. 

Mr. Neison is not content with what he has already said to 
establish the serious imperfection of Delaunay’s method, but 
fortifies himself in the belief of it by a new line of argument 
(pp. 432-437), where he gives his conception of the essential 
nature of Delaunay’s transformations. But his argument is 
fatally vitiated because he will have it that the transformations 
in question are rigorously linear in their operation. Thus, to 
illustrate, suppose Delaunay has 

Operation 1. 

Eeplace a 0 by a 1 -h/ficq, e, &c.) 

Operation 2. 

Eeplace a x by a 2 + f 2 {cc 2 , e,, &c.) 

(I use the subscripts, which Delaunay has not, that my meaning 
may be clear.) According to Mr. Neison’s way of looking at 
things, these two operations are equivalent to 

Eeplace a Q by a 2 -f e 2 , &c.) +f 2 (a 2 , e 2 , &c.) 
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Thus he fails to see that Delaunay intends the a l9 under the 
functional sign / 1? to be eliminated by the substitution of 
Operation 2 , as well as the a l which is outside of it. In conse¬ 
quence he misses all the terms which are of the order of the 
product of f y by/ 2 . 

How, suppose that fy belongs to an operation which is con¬ 
cerned with solar perturbations, and f 2 to one concerned with 
planetary perturbations. Then Mr. Heison, by his erroneous 
interpretation of Delaunay’s processes, fails to get some terms 
of the order of the product of the masses of the Sun and 
planet, which, nevertheless, Delaunay has. How, these are the 
very terms Delaunay is accused of neglecting. And. what is 
sufficiently singular, Mr. Heison appears to regard the symbols 
a, e, &c., which are under the functional signs fy, / 2 , <fcc., as 
having everywhere throughout the whole series of operations the 
same signification, and as being absolute constants ; so that, for 
him, all the fs are explicit functions of the time. 

There is another way in which Mr. Heison’s error may be 
illustrated. Suppose we write one of the differential equations 
of the Moon’s motion in rectangular co-ordinates, thus 


dt 2, dx dx dx dx 


a cm o 

dx 


n cm x 

dx 


&c., 


where 0 denotes the potential of the force exerted by a centro- 
baric Earth 5 and the portion of the disturbing function due to 
solar action has been broken into a number of parts R (0) , e / It (1) , 
6 /2 R (2) , &c., severally proportional to the various powers of the 
solar eccentricity e !; and / 3 R 0 is the portion due to the figure of 
the Earth, (3 being a constant which measures the deviation of 
the Earth from a eentrobaric body ; in fine, m" R 1? is the portion 
due to the action of a planet whose mass is m 7/ . Then 
Delaunay’s way of proceeding is very similar to this : he first 
ascertains what would be the expressions for the Moon’s co¬ 
ordinates were R (0) the complete disturbing function, by making 
variable the a, e, y, Z, g, and ii which appear in the elliptic 
formulae; he then transposes ft® over to the left member of the 
equation, and the potential of the principal force is now no 
longer £3 0 but £3 0 + Jt (0) ; he then proceeds to treat e 7 Tt (1) as if it 
alone constituted the whole of the disturbing function, using the 
elements a , e, y, l, g , and 1 %, which stand in his last expressions 
for the co-ordinates as variables, not those which belong to the 
elliptic expressions. When this is done, e 7 !!® is transferred to 
the left member, and the potential of the principal force is now 
£3 0 -f R/ 0) -f e'R (1) , and the work is continued as before. 

How, Mr. Heison admits the legitimacy of all this as long as 
we are dealing with the portions of the disturbing function 
which arise from solar action; but says that, the moment we 
arrive at the term all changes. Then certain ghosts, as 

it were, of the portions R (0) , e 7 R/ l) , &c., unbidden return to the 
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right member and trouble tbe portion Thus we have 

the strange spectacle of forces figuring at once as principal and 
as disturbing. Mr. Stockwell made a precisely similar objection 
to my elaboration of the inequalities due to the figure of the 
Earth, which was disposed of by Prof. Adams in a single 
sentence. 

If all this be true, what becomes of the assertion, often 
reiterated, that, when the differential equations are written 
down, all the rest is a pure question of analysis P On Mr. 
Heison’s and Mr. StockwelTs view, the analyst, who does the 
integrating, needs an astronomical or mechanical prompter at his 
elbow to inform him of the exact physical import of the 
constants p or m /; , otherwise he will infallibly’' go wrong. 

Wa?lvihifn;i. 

1886, Oct. 14. 


{In Kepler's Problem, By Robert Bryant, B.A., B.Sc. 

The solutions of this problem are usually given with a view of 
obtaining the eccentric anomaly directly from a series expanded 
in powers of the eccentricity, the application of which is in 
general very laborious. 

The use of the differential formula 

d E _ I 

c/M i—e cosE* 

when an approximate value of E in a planetary orbit is known, 
involves less labour, while the results converge with great 
rapidity. 

The two following series are given as readily offering the 
required approximate value of E in a planetary orbit, when one 
application of the differential formula usually gives the eccentric 
anomaly with sufficient accuracy. The series are so simple and 
so readily present themselves that they have probably been given 
before, but I have not met with them. 

We have M=E — e sin E with M and e, given to find E. 
Lagrange’s theorem may be concisely expressed thus : 

If 

y = s+ x<t> (y), tlien/(!/) = 2^ 1 'K 2 )!”/'<»} ’ 

the interpretation of which when n = o is obvious. 

If we expand sin E and cos E, we obtain 

sin E = sin M + e sin M cos M + 6 sinM(2 — 3 sin-M) 4- — sinMcosM(3 — 8 sin 2 M) 

2 3 

+ —ein M (24—136 sinM+ 125 sin 4 M) + .... 

24 

cos E = cos M - e sin 2 M - sin 2 M cos M — 2 ~ sin 2 M (3 — 4 sin 2 M) 

2 3 

t fit 

4-4—sin 2 M cos 11(12 — 4 sin 2 M) + .... 


© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 


Downloaded from http://mnras.oxfordjournals.org/ at Florida International University on June 10, 2015 







